General relation between drift velocity and dispersion 
of a molecular motor 
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We model a processive linear molecular motor as a particle diffusing in a one-dimensional periodic 
lattice with arbitrary transition rates between its sites. We present a relatively simple proof of a 
theorem which states that the ratio of the drift velocity V to the diffusion coefficient D has the 
upper bound 2N/d, where is the number of nodes in an elementary cell and d denotes its length. 
This relation can be used to estimate the minimal value of internal states of the motor and the 
maximal value of the so called Einstein force, which approximately equals to the maximal force 
exerted by a molecular motor. 

PACS numbers: 87.16.Nn. 05.10.Gg, 05.40.-a 



I. INTRODUCTION 

Molecular motors are enzymes mediating intracellular 
transport and cellular motion in living cells. The most 
important examples are kinesins (responsible for trans- 
port of organelles) , myosins (which drive muscle contrac- 
tion) and dyneins (involved in cellular locomotion) Q . 

Based on their function and structure, molecular mo- 
tors can be divided into several families. Here we will 
focus on processive linear motors. A motor is linear if it 
moves along a complementary protein fiber (other types 
of motors include protein pumps, involved in transport 
across membranes, and rotary motors). A motor is pro- 
cessive if it can travel a long distance before it detaches 
from the filament. 

A typical example of a processive linear motor is ki- 
nesin. It is a mechanoenzyme which transduces the chem- 
ical energy of hydrolysis of ATP molecules into mechan- 
ical work, which is then used to move vesicles and or- 
ganelles along microtubules. Each step is powered by 
hydrolysis of one ATP molecule. The main theoretical 
challenge is to explain how the chemical energy of ATP 
is transformed into practically unidirectional motion of 
protein motors. 

Theoretical modeling of the motion of motor proteins 
has been mainly based on two approaches. The first one, 
called a thermal ratchet model, regards a motor protein 
as a Brownian particle diffusing in several periodic, asym- 
metric potentials which are alternatively switched on and 
off in stochastic time intervals [||. In the following we 
employ the second approach, which is based on an as- 
sumption tliat the motion of a motor can be modeled 
as a sequence of transitions between discrete mechano- 
chemical states of the motor, which are represented as 
nodes of a linear lattice ^, |[ |^. A segment of a mi- 
crotubule corresponds to an elementary cell of the lat- 
tice, and the number of sites in an elementary cell, N, 
is equal to the number of different internal states in a 



full mechano-chemical cycle of a motor. For kinesin the 
lattice constant d = 8 nm @, |, |7| and iV > 4 H. 

Any theoretical model must account for experimental 
results. Here we will be interested in an approach which 
attempts to determine internal properties of molecular 
motors based on the statistical analysis of their motion 
as observed in experiments. The fundamental quanti- 
ties employed in this method are the mean velocity V = 
{x{t)} /t and the dispersion D = {{x{t) - {x{t)) f)/2t. 
Having determined these two parameters (in the limit of 
large t), one can find the so called randomness r defined 
simply as [|, |] 



r = 2D/Vd. 



(1) 



This quantity can be used to estimate the number of 
internal states of a molecular motor. The idea is simple: 
if transitions constitute a Poisson unidirectional process 
with exactly the same transition rates fc, then r = l/N 
. This idea was extended by Fisher and Kolomeisky 
^, who made a conjecture that in a general case of 
a linear chain with bidirectional Poissonian transitions 
there is 



N > 1/r, 



(2) 



which enables one to determine the minimal number of 
internal states in a full mechano-chemical cycle of a mo- 
tor. 

Another important quantity which can be directly re- 
lated to V and D is the so called Einstein force 



V 



(3) 



where T is temperature and ks denotes the Boltzmann 
constant. This quantity is a linear-response estimation of 
the force exerted by a molecular motor |^. Conjecture 
(^) can be now written as 



Fe < 2kBT 



N 



(4) 
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The maximal value of the Einstein force is thus a sim- 
ple function of three parameters: temperature T, a step 



2 



length d and the number N of different internal states of 
a motor. 

In a recent paper we have given a rigorous proof of 
relation for linear chains with arbitrary Poisson tran- 
sitions between nearest-neighbor sites. However, that 
proof is rather lengthy and complicated. Here we present 
a different argument justifying and (|^). Not only is 
it much simpler, but it is also more general, as it can 
be used for systems with transitions between arbitrary 
lattice sites (provided the lattice is periodic). 

The paper is organized as follows. In Section || we 
present a mathematical definition of the model. Section 
III contains the proof of equation and Section |^ is 
devoted to the discussion of our results. 



II. MODEL 

We consider a one-dimensional lattice with its sites lo- 
cated at Xn, n (z Z. At time i = we put a particle at 
site xq = and let it jump between the nearest-neighbor 
lattice sites. Transitions are assumed to represent a Pois- 
son process governed by the master equation 



at regular intervals r and are characterized by probabil- 
ities kmnT (kmriT IS the probability of jumping from n 
to m). It has been rigorously shown that the drift 
velocity V and dispersion D calculated in the original, 
continuous- time system are related to the drift velocity 
and dispersion of its discrete-time counterpart 
through general formulae 



V ^v"- 



D = D^ 



(7) 



dP{n, t) 
dt 



= fc+_iP(n-l,t)+fc-+iP(n+l,i)-(fc++fc-)F(n, 

(5) 

where P(n^ t) denotes the probability of finding the par- 
ticle at site Xn at time t and fc^ > are the (constant in 
time) transition rates from a site Xn to Xn±i- We assume 
that the system is periodic in space with a period N > 1 
and a lattice constant d > 0. We do not demand that 
the distances between consecutive lattice sites, — a;„, 
should be all equal to each other. Our goal is to prove 
(S) for any choice of A^, fc+ and k~. 



III. PROOF OF RELATION (§) 



Relation (0) is equivalent to 



D > 



2N' 



(6) 



We will prove this inequality by fixing the values of V , d 
and N and finding the minimal value of D as a function 
of d and transition rates , n = 1, . . . , A. Since V/ d 
and D/d^ depend only on transition rates k^ the 
problem of proving (^ reduces to the one of minimizing a 
rather complicated function of 2N nonnegative variables 

^" ' n 

Probably the simplest proof of inequality (^) is based 

on a very general and nontrivial property of diffusive lat- 
tice systems [|o| . Suppose we have a diffusion process on 
a periodic lattice with Poissonian transitions between its 
sites and let kmn > denote the transition rate from site 
n to m. There exists a time interval r such kmnT < 1 for 
all n and m. We can thus define another stochastic pro- 
cess on this lattice, in which transitions can occur only 



The problem of finding the combination of transition 
rates k^ which minimize the value of dispersion D for 
some fixed values of V and d in the continuous-time pro- 
cess is therefore equivalent to the same problem posed 
for the corresponding discrete-time system. However, 
the latter problem is much simpler. By definition dis- 
persion cannot be negative. Moreover, — for 
two particular choices of transition probabilities fc^r: ei- 
ther when k~T = and k^r = 1 or when k~T = 1 
and fc+T = 0. In these two cases the motion of a diffus- 
ing particle is actually unidirectional and deterministic 
— hence the value of dispersion — 0. Clearly, 
assumes the minimal possible value ordy in these two par- 
ticular "degenerated" cases where the stochastic, bidirec- 
^^ tional process turns into a unidirectional, deterministic 
one. Non-vanishing transition rates in the corresponding 
continuous-time process are equal l/r, and in this case 
V = d/Nr and D = {d/Nf /2T = Vd/2N @. The 
minimal value of dispersion D is thus equal to Vd/2N , 
which completes the proof of (^. 

It is worth noting that not only have we just proved 
inequality (|^) , but but we also showed that V/ D assumes 
its maximal value, 2N/d, if and only if all transitions are 
unidirectional and of the same magnitude, i. e. when for 
E^ny 1 < n, TO < A^ there is 

k+ = k+ > 0, k- = or fc- = k- > 0, k+ = 0. 

It is also interesting to note that in our present ap- 
proach we almost did not employ the assumption that 
only transitions between nearest-neighbor lattice sites are 
allowed. The same line of reasoning could be applied for 
more general systems. The only difference would be that 
in a general case there would be several choices of tran- 
sition rates for which the auxiliary discrete-time process 
is deterministic. Some of them could be characterized by 
different values of A and we would have to plug the small- 
est of them in Eq. (H). However, in the case of molecular 
motors such complicated models do not seem to be ap- 
propriate, as we expect that their behavior is governed 
by the most probable reaction path [Q. 



IV. CONCLUSIONS 

We have rigorously proved that in diffusive processes 
on periodic lattices with Poissonian transitions between 
its nodes the value of diffusion coefhcient D is bounded 
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from above by Vd/2N, where V is the mean drift veloc- 
ity, d is the length of the elementary cell, and N denotes 
the number of sites in the elementary cells. The present 
approach is much simpler than our previous method em- 
ployed in Rcf. 1^. It is also more general, as it can be 
applied to arbitrary lattice systems, while our previous 
approach was applicable only to linear lattices with tran- 
sitions restricted to the nearest-neighbor sites. 

Our result has interesting applications in the theory of 
molecular motors. First of all it implies that by measur- 
ing experimentally the randomness r = 2D/Vd one can 
determine the minimal number N of internal states of a 
motor in its full mechano-chemical cycle. As a matter 
of fact, this approach was already employed in Ref. 
where it was found that for kinesin A'^ > 4. Our result 
gives firmer basis for this kind of argumentation. 

Our result implies also that the maximal value of the 
Einstein force, Fe, is 2kBTN/d, where T is the tempera- 
ture and kg - the Boltzmann constant. We showed that 
this maximal value is attained if and only if all transitions 
are unidirectional and of the same magnitude. 

Although our approach is valid for any value of N and 
an arbitrary choice of transition rates, in the context of 
molecular motors it still refers to a simplified situation. 



First, we assumed that the motion of a motor protein 
can be reduced to diffusion of a Brownian particle on 
a linear chain of lattice nodes, and other topologies de- 
serve at least equal attention |§, ^, 11, Second, we 
assumed that transitions constitute a Poisson process. 
Although this hypothesis is confirmed by recent exper- 
iments on myosin JlSf , it is well known that transition 
rates with other probability density functions may lead 
to completely different upper bounds for dispersion D 
(and hence the upper bound for Fe and the lower bound 
for r may also change) [|[ ^, ||. Nevertheless, the force 
calculated using all these approximations (4.3 pN) agrees 
quite well with the experimental value of the stalling force 
Fg for kinesin (different experimental techniques yielded 
~ 4-8 pN Jill ). Further work is, of course, required to 
clarify the relevance of the above-mentioned problems. 
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